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Ultra-light scalar fields provide an interesting alternative to WIMPS as halo dark matter. In 
this paper we consider the effect of embedding a supermassive black hole within such a halo, and 
estimate the absorption probability and the accretion rate of dark matter onto the black hole. We 
show that the accretion rate would be small over the lifetime of a typical halo, and hence that 
supermassive central black holes can coexist with scalar field halos. 
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I. INTRODUCTION 

The standard assumption concerning galaxy dark mat- 
ter halos is that they are comprised of some weakly- 
interacting massive particle (WIMP). However recently 
there has been increased interest in an alternative pos- 
sibility, that the dark matter halo may be comprised of 
some ultra-light scalar field |, § § |, g § 0, § 1 0, [l] . 
A large number of such scalar particles, all in their 
ground state, can be bound by their self-gravity; the con- 
figurations possess a core radius related to the Compton 
wavelength of the particles in question, and for suitable 
choices of parameters such halos can give a good descrip- 
tion of observed rotation curves [l0| , and optimistically 
may even provide possibilities to alleviate the 'cuspy core' 
H §, [tC], [lj and 'substructure' & § § problems of the 
standard WIMP hypothesis ^ £3^ 

Development of this scenario is at a primitive stage 
compared to the WIMP hypothesis. While it is known 
that the linear theory evolution of perturbations matches 
the standard scenario, and that time-independent equi- 
librium configurations can broadly reproduce desired 
halo properties, the scenario has not been developed in a 
full cosmological setting where halo formation is tracked. 
Nevertheless, what is known so far is sufficiently intrigu- 
ing that the scenario merits further study. 

In this paper, we address one requirement of the scalar- 
field halo model, which is that such a halo must be able 
to survive the existence of a supermassive black hole at 
its centre, as it is widely believed that such black holes 
reside within many or perhaps even all galaxy halos . 
In the WIMP scenario, the angular momentum of the 
individual dark matter particles, combined with their low 
interaction rate, ensures that the capture cross-section 
for halo particles by the central black hole is sufficiently 
small. However, the scalar-field halo regime is markedly 
different; the individual particles do not possess angular 
momentum and indeed are expected to have a Compton 
wavelength upwards of one parsec so that the individual 
particles occupy a considerable volume of space. It is 
important to verify that the halos are able to survive 
the presence of a central black hole if the scenario is to 
remain feasible. 
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This paper is constructed as follows. In Section D, 
we describe the basic steps to model a spherical scalar 
halo without luminous matter. The main intention is 
to provide a simple panorama of the modeling and its 
appealing properties, such as the smooth scalar profiles. 
In Section III, we use two complementary views of the 



interaction of a scalar halo and a black hole: the classical 
Newtonian picture and the semiclassical approximation. 
The latter will give us information about the absorption 
probability and the accretion rate of scalar matter onto 
the black hole, the main result of this paper. Finally, 
we discuss the main results and some points deserving 
further investigation. 



II. SCALAR FIELD HALOS 

We briefly describe a galaxy halo assuming that it is 
made only of scalar field matter. A description includ- 
ingjfor instance, an exponential disk of luminous matter 
[0, pLOfj , would not significantly change the final results. 
Two similar but distinct kinds of scalar field objects have 
been proposed in the literature to explain galaxy halo 
structure: boson 'stars' (comprised of a complex scalar 
field) jy| [T(H [r7| [f^| a nd oscillatons (made from a real 
scalar field) ]fl|, |f9|, KQ|. For simplicity, we restrict our- 
selves to the case of boson stars, though the main results 
can be easily extended to the case of oscillatons. 

The simplest boson stars are those possessing spherical 
symmetry, for which the metric is written in the form 



ds 2 = g^dx^dx" = -Bdt 2 +Adr 2 +r 2 (d6 2 + sin 2 

' (1) 

where A(t, r) and B(t, r) are functions to be determined 
self-consistently from the matter distribution. At the 
classical level, a complex scalar field $ endowed with 
a scalar potential V"(|$|) is described by the energy- 
momentum tensor 



<S>:^. M -g^(<S»°^ a + V)} . (2) 



A self-gravitating boson star is found by solving the cou- 
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pled Einstein-Klein-Gordon (EKG) equations 



dV 
dV_ 



(3) 



where G„ 



is the Einstein tensor corresponding to the 
metric Eq. ([!]), kq = 8irG (we are taking units such that 
c = h = 1) and □ is the covariant d'Alambertian opera- 
tor. 

The EKG equations Eqs. (||) admit solutions of the 
form y/~Ko& — (f>(r)e~ luJt . With such an ansatz, the scalar 
energy-momentum tensor Eq. (|^) and the metric func- 
tions g^u in Eq. ([!]) are time independent. If we now 
search for regular and asymptotically flat solutions, we 
should set the boundary conditions <j)'(r = 0) = 0, A(r = 
0) = 1 and <f)(r = oo) = 0, A(r = oo) = 1, respectively. 
The EKG equations are then reduced to an eigenvalue 
problem; for each central value of the field (f>(r = 0) = <fio, 
it is necessary to determine the (eigen)values of the fun- 
damental frequency w and B(r = 0) = B to find solu- 
tions in which the field has n nodes and satisfy the above 
boundary conditions. 

In principle, we should also impose the boundary con- 
dition B(r ~ oo) — 1. However, the eigenvalue prob- 
lem is further simplified since we can absorb w into the 
metric function B. In this way, u) does not appear ex- 
plicitly in the EKG equations and then it becomes an 
output value determined by uj/m =1/ \jB(r = oo). The 
normalized temporal metric coefficient is calculated via 
g tt = ~{uo/m) 2 B{r). 

According to observations, the gravitational well in 
galaxies is quite weak, which suggests that we should 
seek boson star solutions in the weak-field limit. It is 
then appropriate to choose a quadratic scalar potential 



S I 0- 



This choice is made not 



vm) = 

only for simplicity, as a quadratic potential can also be 
considered an approximation to more complicated ones 
possessing a minimum || @, §j, § |Tl| . 

Using the dimcnsionlcss radial coordinate x — rm, 
the EKG equations become the so-called Schrodinger- 
Newton (SN) equations @, [k], [j]], |lj in the weak- 
field limit: (</>o, —gtt — T <?rr — 1) <C 1- Thus, we need 
only solve the simpler set of ordinary differential equa- 
tions 



(x4>)" = xU4>, 

(xU)" = x<t) 2 , 



(4) 
(5) 



where primes denote derivatives with respect to x. In 
order to clarify the meaning of function U (x) , we take a 
look at the metric coefficients in the weak-field limit, 



g tt si l + U(x)-U ao , 
g rr ~ 1 + xU (x) . 



(6) 



frequency is given by (w/m) 2 = 1 + Uoc, with U(x — 

OO) = Uoo. 

Despite its simplicity, the system above still has to 
be solved numerically, with the different solutions char- 
acterized by, for example, the central value cf>o. As in 
the relativistic case, the solution of Eqs. (^) and (||) is 
an eigenvalue problem; we have to find the one value 
[/(0) = Uq in order to satisfy the boundary conditions 
stated above and to find n-node solutions of the scalar 
field 4>(x). 

To give an order of magnitude estimation of the quan- 
tities involved, the scalar halo models in the literature 
[jl], ||, ^|, fil], [ll| consider an ultra-light boson mass 
m ~ 10 _2 ^V, whose corresponding Compton length is 
Ac = mT x r~j 1 pc. On the other hand, the central ampli- 
tude of the scalar field would be proportional to the gravi- 
tational well in galaxies, and then 0o ~ \Uq\ ~ v 2 ~ 1CP 6 
with v the rotational velocity of luminous matter in 
galaxies (in units of c). 

All information of the properties of the scalar halo is 
contained in Eqs. (0) and (||). Of special interest are 
the asymptotic behaviors of the scalar and gravitational 
fields near the centre. It is easy to show that |2l| 

cf>(x) = cf, Q [1 + (l/6)U Q x 2 + 0(x 4 )] , (7) 
U{x) = Uo + {l/6)4> 2 x 2 +0(x 4 ), (8) 

and so the scalar field remains constant up to radii of the 
order r ~ |CA)| _1 / 2 Ac ~ lkpc. Therefore, the resulting 
self-gravitating object has a smooth central profile up to 
distances much larger than the Compton length of its 
particles. 



III. THE CENTRAL BLACK HOLE 

The geodesies of scalar halos allow massive particles 
to reach the centre of the halo, and in principle the ac- 
cumulation of matter at the centre is not prohibited. 1 
Therefore, a black hole can form in the centre of bosonic 
objects and become a threat to their existence. Such 
a caveat has been recognized before but it is only 
recently that the interaction between black holes and cos- 
mic scalar fields has begun to be studied seriously |24|, . 

Our aim now is to outline the physical consequences of 
the interaction between black holes and the scalar halos 
considered above. For this, we will take the two simplest 
approximations at hand: the classical picture, in which 
the black hole is taken as a central point-like gravita- 
tional source, and the semiclassical picture in which the 
scalar field lives in the curved space-time outside a black 
hole. As we shall see below, these approximations are 



Hence, the usual Newtonian potential is given by Um = 
(l/2)[U(r) — Uoo], while the value of the fundamental 



1 However, it has been shown that rotating Newtonian boson stars 
could provide extra repulsive forces at the centre p^ |, which 
suggests that the inclusion of rotation could avoid the excessive 
accumulation of matter at the centre of scalar objects. 



3 



complementary and can be matched into the scalar halo 
picture of Section 0. 



A. The classical picture 



Taking into account that the Schwarzschild radius 



2GM] 



bh 



9.57 x 10~ 14 pc 



Mi 



bh 



Mr. 



(9) 



of a central black hole is much smaller than any of 
the typical length scales present in realistic scalar halo 
(e.g. the optical radius r opt or the scalar Compton length 
to -1 ; see Ref. |l4| and references therein), we can deal 
with it within the Newtonian regime. 

In the classical picture, the total gravitational poten- 
tial is the superposition Un = (1/2) [[/(r) — Uoo — r s /r], 
where U(r) is the scalar self- gravitational well and (r s /r) 
is the gravitational field of the black hole, which can be 
seen as a solution in vacuum. Thus, we need only modify 
Eq. (J4|) to include the gravitational influence of the black 
hole 



(x(j))" = x (U — mr s /x) 4> , 



(10) 



which resembles the Schrodinger equation in a Coulomb- 
like potential ~ 1/r. We can still construct regular solu- 
tions for the scalar field and the other metric functions, 
but we need to change the boundary condition of the 
radial derivative of the scalar field at r = to 



</>'(0) = -</> mr s /2. 



(11) 



The other boundary conditions remain the same. 

In this classical picture, we notice that the black hole 
only affects the behavior of the field at small r, but the 
scalar profile is still regular. At large distances, the scalar 
profile is unperturbed by the presence of the central black 
hole. In other words, in the Newtonian regime the exis- 
tence of the scalar halo is not threatened by the central 
gravitational source. 



B. The quantum held theory picture 

The approximation we now make is to consider that 
the scalar field lives near the horizon of the black hole in 
a fixed Schwarzschild background 



ds 2 = ~g(r)dt 2 + ^-+r 2 (d6 2 + sin 2 
9{r) 



(12) 



where g(r) = 1 — r s /r, and then its properties are de- 
termined by the field theory in such a curved space- 
time. This is reasonable since, as stated above, the self- 
gravitating effects of the scalar field appear only at dis- 
tances of order r 3> mT x r s . 



Recalling that we are working with a quadratic scalar 
potential, an s-scalar wave 2 obeys the Klein-Gordon 
equation in metric Eq. ( |l2] ) 

1<9 2 $ 
g ~W 



1 d 



dr 



r 9ir- 



dr 



TO $ . 



(13) 



with the corresponding equation for the complex con- 
jugate field $*. Eq. ( p"3| ) is separable in the form 
v / Ko < i>(i, r) — <fi(r)e~ lmt , where we have set u> — to, an- 
ticipating the classical result in which the fundamental 
frequency does not appear explicitly. The change of vari- 
able preserves the notation of section ||. 

At this point, it is convenient to take the Schwarzschild 
factor g(r) itself as the independent variable. Then, the 
differential equation of </>(r) near the horizon is 

5 2 ^+#' + TO 2 r 2 (l-«?)- 3 = O, (14) 

where prime denotes derivative with respect to g. 

The ingoing solution of Eq. ( p"4| ) is given, around g = 
(r — r s ), in the series form (found using the computer 
algebra package maple, www.maplesoft.com) 



$(u,r) = $ (0) (v,r) 



1 



(TO7* S ) 2 ( P n 
n=l 



imr s Q n )g r 



where 



-im[v—r—r s ln(r/r s )] 



(15) 



(16) 



Here v = t + r* is the advanced time coordinate defined 
via the usual Kruskal coordinate r* = r + r s ln(r/r s — 1) 
and we have used the relationship 

g r *( r ) = e r.-r-r s ln(r/r B ) _ ^ 

The coefficients P n , Q n in Eq. ( |l5| ) have the complicated 
form 

3 

Pi 



Pi 



Pa 



1 + 4m 2 r 2 ' 

3 2 + 5m 2 r 2 + 6m 4 r^ 
2 2 (1 + TO 2 r 2 )(l + 4m 2 r 2 ) ' 
1 40+110TO 2 r 2 + 15lTO"M + 36TO 6 r? 



2 2 3 2 (1 + (4/9)TO 2 r s 2 )(l + TO 2 r s 2 )(l 
6 



1 + 4m 2 r 2 ' 
3 2 



2 2 (l + TO 2 r 2 )(l +4m 2 r2)' 
1 80 - 266m 2 r 2 - 319mM 



108m 6 r Q 6 



2 2 3 3 (1 + (4/9)m 2 r 2 )(l + TO 2 r 2 )(l + 4w 2 r 2 ) ' 



2 This is the scalar wave with lowest angular momentum I = 0, 
and hence also the lowest energy. This condition is satisfied for 
the scalar halos considered so far, which are supposed to form 
a (ground state) Bose condensate. The results could be also 
applied for the case of cosmological scalar fields at late times. 
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In the particular case in which mr s <C 1, we can approx- 
imate P n iQ n by their leading terms. Then, we find the 
approximate expressions 



Pn 



2^(n+l)(n- 



n 



(18) 



With these approximate formulas, the sums in Eq. ( O ) 
can be written in terms of known functions, which indi- 
cates that the series diverges for g —> 1 (r — ► oo). 

However, we find that for distances mT 1 > r » r s (for 
which we can neglect (r s /r) 2 and higher-order terms), 
the radial equation for the scalar field becomes 



+ 77i r s — = . 
r 



(19) 



where primes now denote derivatives with respect to r. 
The new solutions arc of the form 



(r) = r- 1 ' 2 [CJi + DYi] [2y^ 



(20) 



where J and Y are the Bessel functions of the first and 
second kind, and C and D are arbitrary constants. 

The overlap region between the two solutions Eqs. jl5| ) 
and (|2C[ ) is mT 1 » r ^> r s . As we said above, using the 
approximate formulas (|l8|), we can estimate the sum of 
the series in Eq. Jig), For example, if r = 10 3 r s , 

$ ~ $ (0) (v, 10 3 r s ) [1 - (mr s ) 2 (512 + Uimr s )] . (21) 

The factor (mr s ) highly suppresses the contribution of 
the series in Eq. (|l5|), so that we can safely approximate 
the radial part of the latter in this region as 



-) ~ 1 — imr 2 /r . 



(22) 



On the other hand, for distances r <C m 1 , Eq. (p 
reduces to 



(r) ~ (m 2 r s ) 1/2 C 



{m 2 r s )- l / 2 D 



(23) 

Notice that we have included a first-order term in 



Eq. (23), just to show that the next-to-order correction 
is simply the Coulomb-like one, which coincides with the 
classical picture above Eq. flll)). 

Matching Eq. (^3|) onto Eq. in the overlap region, 
we find 



D 

C 



= iTr(mr s ) , 



(24) 



which gives the absorption probability of an I = spher- 
ical wave as [E6| 



r = i - 



D in/2 



1 



D . 



-in/2 



47r(mr s ) 3 , 



(25) 



where we have again assumed that mr s <C 1. The inter- 
pretation of r is that it gives the fraction of the ingoing 



wave, and hence the fraction of the incoming particles, 
which is absorbed by the black hole. 

The last result indicates that for typical values mr s ~ 
1CP 7 , we have T ~ 10~ 20 which implies that the absorp- 
tion of the scalar field is negligible and that, from the 
semi-classical point of view too, a central black hole and 
a scalar halo can be put together. Eq. ( p5f ) coincides 
with previous calculations, which also indicate that the 
absorption probability of higher Z-modes is further sup- 
pressed by a factor of the order (mr s ) 21 p7| . 

Summarizing, we can say that the solutions of the 
scalar halo are given by Eqs. © and © for r > m , 
by Eq. ( p0| ) for to -1 > r r s and by Eq. (JTsJ) for r ~ r s , 
with the absorption probability Eq. ( p5| ) calculated in the 
overlap region r s r < mT x . Formally speaking, the 
three different solutions are well matched to each other if 
we multiply Eqs. (|l5|) and (^0|) by the central amplitude 
</>o calculated for the scalar halo in Eqs. ( fio| ) and (||). 
Since the latter is just an overall factor, the absorption 
probability Eq. (25]) remains the same. 

Observe that the second solution in Eq. ( p0[ ) could have 
been obtained within the classical picture in Eq. ( |Io| ) , but 
it was not taken into account because it diverges at the 
origin and our purpose was to construct regular solutions. 
But, as we have seen in this section, this second solution 
contains the information of the interaction between the 
black hole and the scalar field, since it is through it that 
we obtained a non-null absorption probability. 

Another important issue that can be calculated is the 
accretion rate of the scalar field into the black hole by us- 
ing the formula Eq. (3.1) in Ref. p3| . The scalar energy- 
momentum tensor should be written in the new variables 
(y,r), and then we obtain for the flux of Killing energy 
across the horizon 

dM/dt = iirr 2 x T vv (v,r s ) = (2G)" 1 (cb mr s ) 2 , (26) 

in which we have included the overall factor <pQ. Using 
typical numbers (j)Qmr s ~ 10 -13 , the accretion rate is 
quite small, dM/dt ~ 10 _14 Af Q y _1 , a result that is con- 
sistent with the small absorption probability given by 
Eq. (p|). 



IV. CONCLUSIONS 

We have analyzed the impact of a central supermassive 
black hole on galactic halos comprised of ultra-light scalar 
particles. From simple physical grounds, we should ex- 
pect that the accretion rate of a scalar halo onto a black 
hole is small, since the boson particles cannot 'fit' into the 
horizon due to their large Compton length. Besides, the 
absorption probability should be proportional to the ratio 
of the effective 'area' of the two objects, and hence pro- 
por tional to (mr s ) 2 as happens for massless scalar fields 

We found that the absorption probability is decreased 
by an extra factor mr s , which assures the coexistence of 
the bosonic halo and the central black hole. For this, 
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we showed how to construct consistent and regular solu- 
tions on different scales. In addition, the accretion rate 
of scalar matter onto the black hole is so small that the 
matter absorbed by the black hole is much less than a 
solar mass in the whole lifetime of the Universe. On the 
other hand, this result would indicate that the current 
observed accretion in galaxy black holes would not be 
due to matter provided by a scalar halo. 

We have only investigated the equilibrium state of a 
relaxed scalar halo and a central black hole, and it would 
be interesting to have a more dynamical view studying 
the formation (simultaneously or not) of the two objects. 
This would require of the evolution of the full Einstein 
equations, which is well beyond the scope of this paper. 

A related issue is the interaction of primordial black 
holes with cosmic scalar fields, prior to the gravitational 
collapse of density perturbations, as recently outlined in 
Refs. |2j|. For a cosmic scalar held endowed with 



a quadratic potential, the accretion rate would also be 
given by formula Eq. (|2^) , and then the held would have 
the oscillatory behavior Eq. ( |l5| ) near the black hole hori- 
zon. That is, the mass of the bosonic held still prevents a 
strong interaction between black holes and cosmic scalar 
helds. Other kind of scalar potentials would lead to more 
interesting pictures [^5] . 
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